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Abstract 

We consider different models of stochastic dissipative equations and theoretically compute the 
probability distribution functions (actually the associated large deviation functions) of the time 
averaged injected power required to sustain a nontrivial stationary state. We discuss the results 
and in particular draw from our results some general features shared by these distributions in 
realistic dissipative systems. 

PACS: 05.40.-a: Fluctuation phenomena, random processes, noise, and Brownian motion. 
02.50.-r: Probability theory, stochastic processes, and statistics. 

1 Introduction 

Among out of equilibrium statistical systems, strongly dissipative mediums have a very particular 
status. Usually, statistical theories of non equilibrium systems deal with tools inspired or inherited 
from thermal equilibrium, a scheme that Onsagcr theory cxamplifics especially well [1, 2]. That such 
a continuity may be possible in numerous cases is always due to the fact that genuine equilibrium 
keeps a certain relevance, thanks to a time decoupling which induces partial equilibrations [3]. 
Even in recent theories of nonequilibrium statistical mechanics aiming at describing systems as 
complicated as glassy materials, equilibrium concepts are generalized and successfully adapted, 
since again time separation is there at work [4] . 

In contrast, strongly dissipative systems have a full rest natural "fixed point" or equilibrium 
(something like a zero temperature situation). To avoid such a state and explore the richness 
of their dynamics, the system has to be fed continuously with energy, by means of an external 
action which holds it in a non trivial stationary state (if the injection mechanism has itself some 
stationarity properties). As a result, the statistical state (i.e. the stationary measure) reached is 
very far from any concept of (Boltzmann-likc) equilibrium, and new approaches must be followed. 
For instance, granular matter looks more or less like a gas when sufficiently shaked; actually, the 
so-called "granular temperature" defined as the mean kinetic energy per particle, is a rather hazy 
notion, since its values are not unique when the system is not monodispcrse [5, 6], in a clear violation 
of a basic requirement for a temperature- that is, the equilibration of temperatures of any pair of 
subsystems at equilibrium (or cquivalently, the energy equipartition principle). 

Therefore, these systems, sometimes termed "far from equilibrium", do not belong to the tra- 
ditional statistical physics (as far as we are concerned with relevant degrees of freedom: of course, 
the Navier-Stokes equation can be derived from local thermodynamical equilibrium considerations; 
however, a turbulent velocity field obeys a dissipative dynamical equation (NS, say) and its statis- 
tics has nothing to do with standard statistical physics) and consequently requires alternative 
approaches. 

Some years ago, such an original approach was proposed [7], which tried to look at the problem 
of dissipative systems from a global and phenomenological point of view. In particular, it was 
recognized that the common feature of all dissipative systems is the energy evolution equation 



* e-mail: farago@lcp.u-psud.fr. Tel: (33)-l-69-15-41-95 



1 



E = I — D, which explicites and summarizes the physics of these systems from a "macroscopical" 
(or experimental) point of view : two channels of energy flow compete, one corresponding to the 
energy injected by the external operator (J), the other related to an inner dissipation term (D, 
easily recognized, since it is in general a volume term proportional to a dissipation coefficient). 
Of course, the physics is actually more complicated, since these energy flows induce complicated 
structurations of the local fields, which themselves determine in fine the statistical properties of 
I and D; but it was thought that if something common to all dissipative systems could be once 
firmly asserted, this will be done necessarily using that kind of general evolution equation, which 
transcends all peculiar details of each experimental situation. 

Experimental measurements were first performed in a turbulent von Karman flow, and statistical 
properties of I were extracted and analysed in stationary regimes. Surprisingly enough, it was 
shown that the fluctuations of I are non gaussian (despite its a priori "extensive" (on the surface) 
nature), and decay abnormally slowly as Reynolds number increases (i.e. more slowly that a normal 
regression ~ 1/Re 1//2 ). The analysis of these experiments is still the subject of active research and 
debates [8], and the physical interpretation of the observed scaling is up to now not fully completed. 
Similar studies were afterwards made on various systems including granular gases, shell models of 
turbulence, self-organized critical systems, and led to interesting developments: for instance, a very 
convincing way to estimate the number of effective degrees of freedom in a shaked granular gas has 
been proposed [9]. 

From a theoretical point of view, exploring the properties of / and D is not easy, for the systems 
considered are characterized by invariant measures generally unknown. In a preceding paper [f 0] , we 
studied a simple model ( "zero-dimensional" ) of dissipative system, where the statistical properties 
of / and D were calculable. More precisely, we considered the following stochastic equation 

x + jx + V'(x) =ip(t) (1) 

where ift is a white noise and V(x) a potential. We read this equation as the dynamical evolution 
of a single dissipative coordinate (x + r yx + V'(x) = 0) subjected to an external forcing ifj(t) (the 
stochasticity of the forcing mimics the chaoticity devclopped by realistic dissipative systems when 
driven vigorously enough out of equilibrium; this ad hoc choice allows exact computations; we 
verified that purely chaotic dissipative systems (like a periodically forced asynchronous pendulum) 
gives qualitatively the same results, as far as time-averaged observables are concerned (see below); 
however, some interesting studies and comparisons could be performed with deterministic chaotic 
systems) . 

The energy of the coordinate is naturally defined as E = x 2 /2 + V(x), and is conserved in the 
absence of noise and dissipation. This energy obeys the evolution equation 

E = i>x- 7± 2 (2) 

where the above-mentioned structure E = I~D is clearly apparent. The properties of instantaneous 
quantities of such systems are not complicated to compute, but we were interested in calculating 
the distribution of partially time-averaged quantities like 

i r t+T 

e=- dt'I(t') (3) 
T Jt 

Initially, the focus on these time-averaged variables was motivated by the so-called Fluctuation 
Theorem [ff, f2, 13], but beyond this particular debate, this partial time averaging presents a 
major advantage for analysis: it smoothes the short time scale dynamics, and gives the possibility 
to highlight phenomena occuring at slow time scales, broadening the analysis based only upon non 
dynamical considerations. As an example, the systems considered in [10] display a great variety of 
different dynamics, according to the form of the potential V(x): confined, unconfined, activated, 
etc. . . Despite these great qualitative discrepancies, we showed that the probability density function 
of £ is asymptotically the same in the limit of large r (to be precise: lim i logP(e) = /(e) with / 
independent of V; this / is called the large deviation function associated to e). In this result can 
be seen that peculiar details of the short time dynamics of energy have faded away through the 
averaging. We also remarked other interesting facts concerning the pdf of injected and dissipated 
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power (averaged during a time r): the large deviation function of the injected power is curiously 
dependent on the initial conditions set of the system at the beginning of the averaging window: if 
the distribution of e is measured, starting always with the same status for the system (referenced 
here with the initial conditions of position and velocity (xq,vq)), the large deviation function of 
the injected power is the- same as the- dissipated power; conversely, in the permanent regime, that 
is when measurements of e are performed the (statistical) stationarity 1 being reached, rare energy 
fluctuations cause a singularity to arise in the large deviation function of the injected power (but 
not in the dissipation). 

These unexpected properties of the statistics of e rise immediately the question of their gener- 
ality In the models we studied in [10], we chose for ijj(t) a Gaussian white noise (i.e. (ij){t)ij){t')) = 
T5(t — t')), to make the computations as simplest as possible. The purpose of this paper is to study 
generalizations of the systems previously considered, and discuss whether the properties of e yield 
new concepts, possibly useful for more realistic systems. 

In a first part, we recall the results obtained in [10]. Then we present some generalizations 
of the simple model (coloured noise, non linear friction), and discuss the properties of the large 
deviation functions of e associated with them. The details of the calculations are entirely postponed 
in appendices. 



2 Injected and dissipated power in Langevin equation 

In [10], we studied the system 

x + jx + V'(x) = ip(t) (4) 
&(t)i>(t')} = 2r6(t-t') (5) 

which was thought as an intrinsic dissipative system (driven by x + 72: + V'(x) = 0) driven away 
from equilibrium by a Gaussian white noise. Of course, it is not the common interpretation of 
a Langevin equation, which describes usually the thermalization of a particle in a fluid; in that 
context 7± and tp are two different faces of the same action. There we abandoned that reference, 
and in particular made no citation of the Einstein relation T = "/ksT. We studied the statistical 
properties of 

£ = - [ dt'x{t')^{t') (6) 
T Jo 

the injected power averaged over a finite time interval of length r. This task is a priori complicated, 
for e involves explicitely the dynamics. Nevertheless we succeeded and computed the distribution 
Prob(e) for any value of r and any potential (in fact any potential non repulsive at 00). The 
distributions are of course potential dependent, but an prominent feature emerged in the limit of 
large r. We know that at large r, the distribution obeys the large deviation theorem: 

3/, Urn -logP(e) = /(e) (7) 

r — >oo T 

/ is called the large deviation function associated with e. This property is sometimes uncorrectly 
noted P(e) ~ exp(r/(£)). We demonstrated that this large deviation function is in fact independent 
of the potential V(x), what is rather unexpected, since the underlying dynamics is on the contrary 
very sensitive to the form of the potential. We found that if the distribution of the initial energy 
(i.e. at t = 0) is bounded, the large deviation function has the expression 

/(£) = -J(e-l) 2 where e = e/T (8) 

(the maximum of this function corresponds to (s) and is always T in this model since (/) = (D) = 
7(i 2 )). Interestingly enough, this expression is no longer valid if the distribution of initial energy is 

1 characterized by time-independent probability distributions 
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unbounded, which is the case in the permanent regime since Prob{v 2 ) oc exp(-lV/ 2 7)/M- In that 
case, a phenomenon analogous to a phase transition occurs and (8) is valid only above e c = 1/3. 
Below that value, the curve is replaced by a portion of a straight line: 

1 -7(1 -2e) if e ^ 1/3 

(A physical explanation of this singularity is given in [10]). Finally wc computed the large deviation 
associated with the dissipation S = i J Q dt^x 2 (t) and found that it is equal to (8), S = 5/T replacing 
e. 



3 Free and harmonically bounded particle driven by coloured 
noise 

We now consider a slightly different system. The dynamics is given by 

x + 7± + u 2 x = ip(t) (10) 

T / \t-t'\ 
(V>(t)V(i')> = -exp -J l - 

It resembles the previous case (with V(x) = u> 2 x 2 /2), except that the external force acting on the 
system is no longer a Gaussian white noise, but a coloured noise, an Orstein-Uhlenbeck process. 
The first merit of this system is that henceforth, there is no possible confusion in the interpretation 
of our system: the situation described by (10) does not correspond anyhow to a thcrmalization 
process. Actually, a Langcvin equation describing the thcrmalization of a particle in a thermal 
bath with finite correlation time (exponentially correlated) , would have involved a frictional term 

r* 

,-(t-u)/T 



du(T/k B TT )e-^ u > /To v(u) (11) 

Jo 

instead of our single jv [1]. 

The second merit of the system is to provide us with one dimensionless quantity a = 7T0 (in 
the free case co = 0), which is a generalization of the previous case where only one typical energy 
can be constructed. It allows to test the reaction of a given system under sollicitations which differ 
only by a characteristic time. 

The details of calculations are fully postponed in the first appendix. Here we discuss only the 
results obtained. The large deviation function /(e) has again gcncrically two expressions depending 
on the general shape of the initial velocity distribution. If this distribution is bounded at t = 0, 
the large deviation function has the expression (with the definitions (3 = lotq and e = e/r)) 

= 77T- f- £ (« 2 + 1 - 2 /3 2 )^(e) - 6 £ (/? 2 + %)9{e) + (a + 1)(8 + -(a - l)(a 2 - 1 - 4/? 2 )) 
Ioto \ a a e a 

(12) 

where 0(e) is the largest real root of 6 3 - (a 2 + 1 - 2(3 2 )B - 2((3 2 + a/e) = 
An alternative and useful expression of this function is also 

/(e) ~ g(X) + Ac with A such that g'(A) = — e 

g(X) = (a + 1 — ^(A)) with 9(X) the largest positive root of (13) 

2t V / 

(6 2 - a 2 - 1 + 2(3 2 ) 2 - 8f3 2 9 - 4(/3 4 + a 2 - 2/3 2 + 4rA 7 T 2 ) = 



(remark that g is the Legendre transformation of /). If lj = 0, the previous expressions are 
simplified; for example we get 



gp=o(X) = (a + 1 - V /a 2 + l + 2a v /l + 4rA/ 7 ^ /(2r ), 



(14) 
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and / could be as well explicitely computed, since the inversion of g(X) involves roots of a third 
degree polynomial. Incidentally it can be checked that (13) gives the correct limit for To — * 0: 
owing to the fact that 9 = 1 + 0(tq), the leading order yields 

(6 2 - l) 2 = 8/3 2 + 4(a 2 - 2/3 2 + 4rA 7 r 2 ) (15) 



whence <JV 0= o(A) = ^(1 — + 4T\/~f) whatever the value of j3, as established in [10]. 

The typical shape of that function is given on the figure 1 and is typically asymmetric. The 
principal characteristics of these functions are the location and the curvature of the maximum, and 
the nature of the divergence at e — > and e — > oo. The maximum of / is obtained for the average 
injected power e = e (and correspondingly A = A = since = [</(A) + e]d\/de + A), 

(16) 



+ 1 + /3 2 



It is interesting to note that the presence of a confining harmonic potential is always a hindrance 
to injection of energy into the system. Whether this fact is absolutely general whatever the form 
of the (confining) potential is still an open question. 

The curvature of / near the maximum e can also be explicitely computed, for it is easy to show 
that /"(e) = -l/g"(0). We obtain 

j-w.,^, r 2 a 2 + 3a + 1 + /3 2 

This quantity is relevant as one easily measured (numerically or experimentally) besides the mean 
value. Incidentally, we can form with e and aj a characteristic energy 

_aj_ g»(Q) _T f a-p\l + 2a)-(P \ 
cvrv - e ~ 2 5 '(0)~ 7 V (a + l + /3 2 ) 2 ) ^> 

The form of ^T curv (a, (3) is plotted on figure 2. We shall sec that this characteristic energy T curv 
is in a certain sense an invariant of the energy flow process (see later). 

Let us look now at asymptotic branches. The e — > branch is characterized by a e -1 / 3 diver- 
gence: 



This exponent —1/3 is a novelty with respect to the "pure" Langevin case, where the exponent was 
— 1. This is not very intuitive, but follows closely the decreasing of e with a: it is probably due 
to the fact that / = t/jv is bounded as soon as a ^ 0, which was not the case with a white noise; 
this important change in the shape of f(e) shows that the white noise limit is a singular case, even 
for the function f(e) which is a priori related to a time-integrated observable. Besides, it shows 
that the short time dynamical details continue to play a role in the large deviation function. . . We 
can just remark that this vicinity of e — becomes independent of (3: we can thus conjecture 
that this exponent is potential independent and is controlled only by the injection and dissipation 
mechanisms. 

The other asymptotics corresponds to e — ► oo, and is much more softer than the e = one, that is 
a simple linear behaviour, witness of a singularity of g(X) occuring for a finite negative value A = A . 
The slope of the line is homogeneous to (energy) -1 . In the pure Langevin case, f(e) = —^T((e/T) — 
l) 2 /4e and the slope is simply — 7/4I\ We are thus naturally led to define a characteristic energy 
T s i pe by asserting that the slope of the oo divergence is given by — l/4T s ; ope .From (12), 

T slope = -4« 2 ^ (-(a 2 + 1 - 2/3 2 )C 2 - 6C/3 2 + (a 2 - l)(a 2 - 1 - 4/? 2 ))" 1 (20) 
where C largest solution of C 3 - (a 2 + 1 - 2/3 2 )C - 2/3 2 = (21) 
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The figure 3 gives a plot of ^T s ; 

ope & 

function of a and j3. The behaviour of T curv and T s i ope 
is qualitatively quite similar, what is not so surprising: the typical form of g(X) is shown on figure 
4. If we try for g (in the vicinity of A = 0) the Ansatz g(X) w g(Xo) — C'\A — Ao, we deduce 
that would this Ansatz be exact, we would have T curv — T s i ope . As this Ansatz is not so bad (cf. 
fig. 4), we deduce T curv r* T s i ope . Thus, the large deviation function /(e) has its right asymptotic 
branch slightly constraint by the vicinity of its maximum, but this relationship remains however 
relatively weak (the relative difference does not exceed 45%). Physically, the energy T curv is more 
interesting, since it is constructed with the top of the curve, i.e. with quantities which are easy to 
measure/compute, which is not the case with the asymptotic tails, associated to rare events. 



3.1 Singularity of / in the permanent regime 

From the preceding paragraph, it is clear that unlike the white noise case, the large deviation 
function of e is sensitive to the short time dynamics and to the presence of a potential. Thus, 
the irrelevance of V(x) demonstrated in [10] is specific to the white noise case and it is probably 
hopeless to seek in general for such simplifications. Nevertheless, something already observed in 
[10] is still present here and is probably extremely general: in the permanent regime -that is, 
when the initial conditions are not fixed, but sampled from the stationary distribution-, the shape 
is altered and a negative tail appears below a certain positive value s c . This tail is simply the 
straight line g(X c ) + X c e where A c corresponds to e c . This effect corresponds to the fact that for 
small values of e, the probability of these rare events is no longer dominated by the inner dynamics 
(inside the time interval [0, r]) but by rare and very energetic initial conditions (cf. [10] for details). 
This situation is due to the fact that I and D are intimately coupled to the energy which has a 
conservative character. As a result, this phenomenon must also occur in the models considered 
here, and more generally in any case where the stationary distribution of energy is unbounded. 
As an example, consider the case [3 = (for the sake of simplicity). The computation of the so- 
called "fluctuation term" in (e~ Are ) provides a prcfactor to cxp(— rg(X)) which displays a cut for 
A > A c = 2-f(2a 2 + 3a + 1)/T. The corresponding value for e is 



y/9 + 16a 2 + 24a ( a 2 + i + 2aV9 + 16a 2 + 24a) 3/2 

and tends to 1/3 when tq — ► 0, as expected from [10] This singularity has an interesting physical 
meaning: it is intimately related to the three observables (E, I, D) and their mutual dependence 
and balance (the stationary distribution is a property of the dynamics). The disadvantage of that 
phenomenon is that it is located in a region of quite rare events, and corresponds to a second order 
singularity: as a result, it is not easily observable. One way to make the singularity observable 
could be to consider "very hot" and externally controlled initial conditions, where large initial 
energies be likely; this will certainly shift the location of the singularity. For instance, in the 
pure Langevin case a = [3 = it can be shown that it shifts the singularity from e/T = 1/3 to 

e/r = 1/ \/3 + « 0.414 (in the limit of flat initial condition). Thus, this shift would be anyway 
limited, but could be however sufficient to make the singularity measurable. In that case, precursors 
of this singularity could be observable at finite r by plotting \d 2 (log P)/de 2 : this function must 
display a steep jump at the singularity location. 



3.2 Characteristic time of the energy flow 

Above was defined a quantity related to the curvature at the top of f(e) as a 2 = l/|2/"(e)|. In 
fact a 2 is also given by the following expression: 

dt[(I(t)I(0)) st - (I) 2 t ] (23) 

where (. . .) s t denotes the averaging in the (out-of-cquilibrium) stationary regime. To demonstrate 
that, let us compute (e 2 ) — (e) 2 in the large r limit. We have 

P(e)=p(s).cMrf(s))-p((e)).cMrf"((e))(e (e)) 2 /2] (24) 







where p(e) is the preexponential factor of the distribution. It is readily obtained (e ) — (e) ~ 
l/|r/"((e))|. Expanding the second moment, the result is established. 

In all the stochastic models we considered so far, we always noticed that large deviations func- 
tions associated with energy injection and energy dissipation were intimately related, and even 
equal in the vicinity of the maximum. It is easy and instructive to demonstrate this property 
(which is always valid) directly from the expression of a 2 (we recall that the stationarity implies 
(I) st = (D) st ). Indeed, we have 



dt 



(m)+D(t))(E(0) + D(0))) st -(D) 2 st 



((E(oo)-E(0))(E(0)+D(0))) st 



(E)(D) - (E(0)D(Q)) st + / dt 

Jo 



dt 



(D(t)(E(0) + D(0))) st -(D 



(D(0)(E(-t) + D(t))) st -(D) 2 st 



di[(£>(t)£>(0)) ( 



(D) 



(25) 
(26) 
(27) 
(28) 



where we made use only of time translation invariance, assumed in a stationary regime. This 
property allows us to define some quantities naturally associated with the energy dynamics. We 
already saw T curv which thus can be computed from the fluctuations of the dissipation term, but 
the very physical interpretation of this energy is not clear up to now. Similarly, we define a 
characteristic time r e associated with the energy dynamics, 



D 



(1)1 



(29) 



which presents also the remarkable property to be symmetric with respect to the injection and 
the dissipation. These considerations are extremely general and make use only of the E = I — 
D structure and the time translation invariance. Thus the equation (29) is always true in the 
stationary regime, and r e is a correlation time attached to the full energy flow process. Usually, 
we define the correlation time of an observable as tx = [/ °° dt(X(t)X(Q) — (^) 2 )] / (X 2 — (X) 2 ) . 
Thus, r e is more or less a bare time whence tj and td are constructed as 



Tl 



(I? 



(i 2 ) {iy< 



(30) 



(and a similar one with D). These relations are interesting, since if we imagine that r e is fixed (or 
takes it as the natural time unit), we see that at the injection level, the fluctuations level is enslaved 
by the characteristic time of the injection mechanism. At the other side, this balance exists with 
a dissipation time which is possibly quite different from the former; in that case the fluctuations 
are again constrained to adjust correspondingly. This simple reasoning of course does not take into 
account the fact that the "reference" time r e is itself a property of the established regime. 



As an example, let us look again at our toy-model. Noticing that r e = T curv /e, we get 

1 a 2 + 3a + 1 + (3 2 
Te ~ 7 a + 1 + /? 2 



(31) 



This result is independent of T the strength of the excitation, a fact obvious in our case, since no 
dimensionless number can be made with T. But for more complicated models a T dependence is 
to be expected a priori. We notice also the effect of the confining potential: the oscillation time 
lu only prevents the correlation time r e to diverge when tq is large. This is obviously due to the 
unavoidable oscillation at frequencies near u> in the response function of x. 



4 Anharmonic potentials 

Unlike the pure Langevin model, the addition of an anharmonic confining potential makes the 
computation of f(e) untractable. Actually the mapping to an effective dynamics (as exposed 
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in appendix) is no longer feasible, due to a proliferation of new types of terms in the action. 
However, the previous analysis will be qualitatively the same, and the pdf of e has typically the 
same asymmetric shape, with a right tail driven by a term exp(— Te/T s i ope ) (this term is in general 
multiplied by a prefactor ~ e~ v where v is potential dependent [10]): figure 5(a) shows the quantity 
- logP(e) for a model a = /3 = 1, T = 1 and a double well potential V(x) — 0.25 * x 2 — x 4 . We see 
that the convergence of log(P)/r to f(e) is quite slow, due to correction terms of order log(r)/r. 
The other part of the figure shows the comparison of two pdf for the same value of (r, a), but for 
two different potentials (harmonic (3=1 and double- well) . It clearly shows a strong dependence of 
log(P) (and therefore /) with respect to the potential. The positive tail is less pronounced in the 
nonlinear case, which is a bit surprising, since the parameter was chosen such that the curvature in 
the bottom of the double well is identical to the harmonic curvature. Naively, as in a double well (p 4 
potential the frequency decreases with the energy we would expect a scenario similar to an effective 
harmonic potential with a lower well frequency, and thus, an enlargement of the rare positive 
events. This is not observed, thus demonstrating that probably rare energetic events correspond 
to scenarios where the steep external branches of the potential are often explored by the particle - 
very energetic oscillations have in fine an increasing oscillation. Finally, a second order singularity 
in / is also expected in the permanent regime, although difficult to observe, as easily seen in figure 
5 (a): for r = 50 (whereas the typical dynamical time are of order 1), the difference - log(P) — / 
is not negligible (cf the maximum of the curve), and several millions of statistic steps are however 
unable to sample the vicinity of e = 0. . . To conclude, let us mention that a variational approach 
could be performed on these nonlinear models: using (e x ) ^ e^ x ', we could compute for each A 
the best harmonic potential describing the dynamics. 



5 Nonlinear friction 



There is another way to generalize the original Langcvin equation, namely to change the dissipation 
term 77J to a nonlinear term tp(v), where <p is an odd function, positive for positive v. It is worth 
looking at these models, since we can that way test the influence of the dissipation term on the 
shape of P(e). As an example, let us consider the simplest case 



v + <p(v) = ip(t) with (i/>(t)ip{t')) = 2T5(t - t') 



(32) 



It is shown in the appendix that the function 17(A) can be expressed in terms of the lowest energy 
of a Schrodinger equation: 



g(X)/Ga = - Min ( ( 



whence we get 



m/r = 



Max 

(C\v<p\()=er 

<C10=i 



(33) 



(34) 



Thus, the large deviation function comes in this case from an extremalization principle, but this fact 
is specific to the situation considered. The inspection of (33) is instructive, for we see immediately 
that if <p(v) diverges faster than v, g(X) has no longer a singularity at finite A. As a result, the right 
asymptotic tail of f(e) is no longer a straight line: the precise form of this asymptotics depends 
on the dissipation efficiency And we find obviously that if ip(v)/v — > 00, the right tail of /(e) goes 
faster to —00. To be more quantitative, if tp(v) ~ v v for large v (with v > 1), we can show by 
considering the effective potential in (33) that f(e) ~ — e 2u /( u+1 \ Does it affect the left asymptotics 
as well ? To answer this question, the A — > +00 limit must be investigated. The effective potential 
is deeply changed by this change of sign and is now single well and is equivalent to Xv(p(v). As a 
result we expect a minimum energy scaling like A 1 / 2 [14] (provided <p(v) ~ v near zero). This leads 
to a divergence e" 1 near zero, whatever the precise form of (p(v). However the condition ip oc v 
near v = is crucial, otherwise the left tail of / is affected: if for instance (p ~ v 3 , the zero point 
energy of the Schrodinger potential scales like A 1//3 and we find /(e)~e~ 1 / 2 . 
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This example is quite interesting, since it shows that the dissipation mechanism can affect deeply 
the general shape of /(e), and particularly the form of the asymptotic tails. In particular, the 
behaviour of tp(v) for small velocities dictates the form of /(e) near e = 0, whereas correspondingly 
the large v behaviour of ip is responsible for the large e value. Thus, for realistic situations, a study 
of the asymptotic tails of the e distribution could provide an estimation of an effective (in average) 
dissipation/energy relation: in a turbulent flow for instance, the greater the energy present in the 
system, the greater the dissipation; we argue that the implicit relation between these observablcs 
is partly encoded in the asymptotic tails of /(e). 

6 Conclusion 

In this paper, we studied the injected power distributions of several stochastic models, in order to 
precise the results obtained in [10] and test their possible generality. This study showed that, in 
general, the large deviation functions remain sensitive to the details of the short-time dynamics 
of the system, despite the time averaging: this observation makes the Langevin-likc models very 
peculiar, since in these cases the large deviation function is insensitive to the presence of a pinning 
potential. However, the common scenario of dissipative systems in a permanent regime makes all 
these to share some common characteristics: first, they are strongly asymmetric, which reflects 
the fact that we follow a positive conserved quantity, with a systematic dissipation; secondly, they 
display a second order singularity in the permanent regime which constructs a negative tail in 
the large deviation function associated to the injected power (this singularity is not present in 
the dissipation, since the dissipation is always positive). Third, the large deviation functions of 
injected and dissipated power have the same curvature at their maximum, due to the structure 
of the energy equation. We can even conjecture that in general the two large deviation functions 
are strictly equal in a vicinity of the maximum (it is easy to show it for the models presented 
here, but probably often true). From this equality a characteristic time (or a characteristic energy) 
naturally associated to the energy flow can be defined. It would be interesting to measure that 
time in an extended dissipative system, when the system is chaotic and generates itself its disorder 
and the noisy character of the energy injection (instead of being imposed by the operator): the 
typical time scale is in that case controlled by the dynamics of the system itself (for instance, 
what is the behaviour of T e near an instability threshold ?) . More generally, it would be worth 
now performing similar analysis to extended systems where non equilibrium stationary states can 
produce structurations of the system. 

7 Appendix A: large deviation function for the non marko- 
vian system 

We derive in this appendix the main results of section 3. We consider the injected power averaged 
over a time interval r : 



for the model (10). As r is finite, this observable remains a fluctuating quantity, and we examine its 
probability distribution function 7r(e). More precisely, we are interested in the computation of the 
so-called large deviation function, a function which emerges from the consideration of large values 
of r. It can be proved that for large values of r the probability of e verifies 



(this equivalence is sometimes boldly noted 7r(e) ~ e T ^ E \ i understanding the prefactor). The 
function / is nothing but the large deviation function associated with e. We concentrate ourself 
on this sole quantity, since the prefactor is quite more involved to evaluate, and is physically much 
less interesting. 




(35) 



log7r(e) 



rf(s) 



(36) 
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For convenience, we are going to compute two different types of pdf, namely 7r Wo (e) the proba- 
bility of e, knowing that the initial conditions of the process are Wq = (x(0) = xq, v(0) = Vq, ip(0) = 
4>o), and 7r(e) the probability of e in the permanent (stationary) regime. The latter is simply related 
to the former through 

7r(e) = J dx dv dip Pst(xo,vo,tpo)n vro (e) (37) 

where P st is the stationary distribution of the (correlated) variables (x,v,ip). As ir and 7r Wo arc 
different, their large deviation functions are termed / Wo and / respectively henceforth. 



7.1 Path integral representation of the characteristic function 

The characteristic function of 7r Wo (e) is defined by 

(e- A -) wo = Jde 7r W0 ( £ )cxp(-Ar £ ) (38) 

As soon as this function is known, the original pdf can be retrieved via an inverse Fourier transform 
(in the complex A variable). At the level of the large deviation function, equation (38) gives a rapid 
answer: there exists a function g Wo (X) such that 

(e- AT£ ) W0 ~ cxp[r 9wo (A)] (39) 

x wo r — too 

(in the bold sense !), and this function is given by the Legendre transform of / Wo : 

.9w (A) = / W0 (e) - Ae (40) 
/w (e) = A (41) 

(the inversion is simply: / Wo = g Wo + Ae, <? Wo (A) = — e). All that is a priori correct, but for a point: 
it is possible that the prefactor of e rg in the characteristic function diverges for a particular value 
of A, and that therefore the function be not defined for certain values of A, whereas the function g 
remains defined in this range. In that case, the rapid Legendre inversion above mentioned must be 
replaced by a careful analysis of the integral (38) (see [10] and below). 

Let us give now a path integral representation of ( e ~ AT£ ) w ■ It is quite easy to do that in 
our case, since the noise ?p is nothing but a Ornstein-Uhlcnbcck, that is a process described by an 
ordinary Langevin equation: 

^ + t -V = C(*) (42) 

(C(f)C(f')> = 2^(i-0 (43) 
r o 

where £ is a Gaussian white noise. The path-integral representation of the propagator of this 
process is known [10, 15, 16, 17] : 

P^i,r|Vo,0) = e^x / ptyjexpl--— j / dtU + T^) (44) 

J^(o)=tp \ 41 r Jo \ I J 

As a result, the statistical weight of a particular occurence of the noise, knowing that it originates 
from ipo at t = is given by 

x exp ^— — !jq J dt (r ip + ipj J (45) 
Therefore, we get the probability of a particular occurence of the process x (beginning at wo) as 



Prob 



[z(t)],te [o,< 



X(0) = Xq, v(0) = Vq, V(0) = ljj - JV ~ V (x ) 



-(!+« 



P (~ 4T J dt + + vV "( x )) + (v + l v + V'{x))f^ 
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where we defined a = 7T0, and V{x) = cu 2 x 2 /2. It must be noted that the jacobian of the 
transformation ij) — ► x is equal to exp(7r/2) and is taken into account in the preceding expression. 

From now on, it is easy to deduce an expression for (e _AT£ ) Wo > since the distribution of individual 
paths is at hand: 



\ / Wq 



/ [Dx] exp I 






/ w(0)= wq 



dt [t (v + ■yv + vV") + (v + + V')} — A dt v(v + yv + V') 



(47) 



(with obviously v(t) = x(t)). It is convenient to express explicitcly the final values of ip and (x, v): 



w(0)— wo 



4r 



[Vv] exp -— / dt [t (v + jv + vV"(x)) + (v + jv + V'{x))] 2 — A7 / dt v 2 ) (48) 



where E = ^u 2 + V(x). The essential point of the derivation is now the possibility to find four 
constants U, 7, To, a), such that the difference 

A = / dt ([t (v + jv + lu 2 v) + {v + jv + lo 2 x)] 2 + 4rA7« 2 ) 
Jo 

-U f dt [t (v + jv + uj 2 v) + (v + 7« + u) 2 x)] 2 (49) 



implies only boundary terms. If we choose (U, 7, ro,2) fulfilling (we define of course a = T07, and 

f3 = T Q UJ,P = T Q) 



r 2 = U? 2 



a 2 + 1 - 2f3 2 = U{a 2 + 1 - 2/3 2 ) 
+ a 2 - 2/3 2 ) + 4fA 7 = + a 2 - 2(3 2 ) 



-4 ' 



(50) 
(51) 

(52) 
(53) 



it completes the goal and we can write 



e 2^( 1 + «)-K B l 2 - S 0)-# 



w(t)=wi / jj 

[Dx] exp ( -— / dt [? (v + jv + Q 2 v) + (v + jv + Q 2 x)] * ) (54) 



w(0)— wo 



The very reason we made this transformation is that the remaining path integral is closely related 
to the propagator of a "renormalized" process of type (10), where (r, 70,7,0;) — > (T/U, r , 7, cD): 



w(r)— wi 



(0)=w 



u 



4r Jo 



[2?cc]exp ( I dt [? (v + yv + w 2 v) + (v + yv + w 2 x)Y ) = e ^ (1+Q ° x P(w lt t|w , 0) 



(55) 



As a result, 

< e_A ">w „= «P 



1 + a 1 + 2 



2r n 



2fn 



y (iwi e"^ (B i" B » 2) "^ x P(wi,r|w o ,0) (56) 
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The remaining integral has a nice behaviour for large r, since P(wi , t|wo, 0) tends to its equilibrium 
value P s t(wi). Thus, we get the veritable equivalence 

<e-^> wo ~ exp " ^]) >< / dwi e"*^"^-* x P st ( Wl ) (57) 

whence we extract the function (/ Wo 

1 fl + a l + 5\ , , 

5w„(A) = - — (58) 



2 V. 2r 2t , 

It is worth noticing that the subscript wo is useless, for the function g Wo is independent of wo; we 
abandon henceforth it (for g) in the following. 

Let us define 9 = (1 + S)to/to. We can verify that 9 is the largest positive root of 

(9 2 - a 2 - 1 + 2/3 2 ) 2 - 8(3 2 9 - 4(/3 4 + a 2 - 2/3 2 + iTXjr 2 ) = (59) 

when it exists; if not, (e~ XTe ) diverges. It proves the result (13). 

Similarly, We get the characteristic function of the process in the permanent regime by simply 
integrating over w : 

(e- Are ) ~ exp(r 3 (A)) x /"dw^wo eS^^S x P st (wi)P st (w ) (60) 



(Note that the stationary distributions P are not the same for Wo and Wi, since P refers to the 
renormalized process). We have to remark that the same function g is also shared by the stationary 
process. But it does not imply that the associated / Wo and / functions will be the same. 

7.2 Analytical properties of characteristic functions 

To derive the large deviation functions /, it is required to look the analytical properties of the 
prefactor, since the latter can lessen the analyticity domain of the characteristic function with 
respect to that implied by the sole inspection of g. 

A priori, the function g is analytical over the whole complex A space, except on a cut located on 
(but in general not equal to) R - . This cut begins when the root 9 disappears. In this paragraph, 
we considers only the free case u> = 0, for the sake of simplicity. But the discussion remains valid 
for the bounded case as well. 

Let us first consider the non stationary situation: here, for a fixed value of Wo, the prefactor is 
composed by functions as regular as g itself {g has a cut on A < A_ = — 7/4T) times an integral 
over wj. Thus, limitations to analyticity could arise if for some values of A, this integral no longer 
converges. A lengthy computation gives this integral (when converging) proportional to 



([y/l + a 2 + 2ar\ + a} 2 - l)\J \J\ + a 2 + 2ai] - a + 1) 



(61) 



(with r/ = v/T-i- 4rA/7) which again is as regular as g itself. It is concluded that for the non 
stationary case, the prefactor cannot hinder the Lcgendrc inversion, and that / Wo is really given 
by the inverse Legendre of g. 

Quite different is the stationary case, since now two integrals have to converge. Of that over 
Wi we already proved the analyticity over C \ That implying Wo is completely different, and 
another lengthy computation gives it proportional to 

( v 7 ! + a 2 + 2ar) + a + l)" 1 ((77 + l)(3o + 1 - ^1 + a 2 + 2cm?)) (62) 

This time, a new cut appears, since for A > A c = 2(2a 2 + 3a + l)j/T, the preceding expression is 
not defined, which expresses the fact that the integral over Wo is not defined. This extra cut has 
deep consequences on the large deviation function /. 
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7.3 Large deviation functions 

For the nonstationary case, the usual rule applies, that is / is the Legendre transform of g: 

/w (e) = g{X) + Xs 
9 (A) = -e 

For the stationary case, this expression is not always valid: (63) must give A < A c otherwise there 
is a problem. In that case, the inverse Laplace transform A — > e must be lead carefully [10] but 
the result is quite simple: if no allowed A fulfills (63), the large deviation function for the injection 
considered is simply /(e) = g(X c ) + A c e, therefore a portion of a straight line. Thus, we have 

/(e) = if £ > e c . g4 . 

/(e) = g(X c ) + X c e if e < e c 

and £ c defined as g'(X c ) = —£ c - 

7.4 Nonlinear potentials 

For a general potential V(x), the expression (48) remains true. If V is not harmonic, the renor- 
malization procedure is no more possible, since the action involves, besides boundary terms, the 
terms 

(t 2 v 2 - v 2 (t 2 V" -a 2 -l) + V' 2 + v 2 {t 2 V" 2 + 7 2 - V") - ^-v 3 V" + ^V< 

Thus, there is a "proliferation" of new terms which cannot be balanced by a simple redefinition of 
the constants. In fact, it seems that the integrability of the problem is restricted to the linear case. 

8 Appendix B: nonlinear friction 

Here we sketch the computations pertaining to the section 5. From the dynamical equation (32), 
(e~ XTe } = JlThflexpf-^ j\v + ip(v)) 2 -X J\v + <p(v))v\ (66) 



v, care must be taken of the fact that the Jacobian of the transformation is \ /J" (p'(v)dt. Then 

(e- XTE ) Va = J dv x ^ {vl ' vl) [ T q _ V"] ^(-^^ [(v + v{v)) 2 -2Tv'{v)+WXv V {v)}} 

(67) 

The idea is the same as for the linear friction cases: we seek for a function £(v) and a constant Q 
such that for all v 

tp{v) 2 - 2T(p'{v) + 4T\v<p(v) = £ 2 (v) - 2T£'(v) + Q (68) 

with the "initial condition" £(0) = 0. This is a Riccatti equation, and the standard transformation 
£(v) = —2r£'(v)/£(v) leads to a Schrodinger equation: 

-4r 2 C" + [if 2 - 2T(p' + 4TXvtp]C - QC = (69) 

Owing to the fact that the potential ip 2 — 2T(p' + 4TXv<p is even, that £ never vanishes (since £ must 
be defined Vt>), we deduce that Q is necessarily the lowest eigenvalue of the Schrodinger problem 
and C the associated eigenfunction. 
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Following the trick of the previous Appendix, we easily deduce that g(X) = —Q(\)/4T, that is, 
using variational properties of Schrodinger eigenvalues: 



To get /, we must take the Legendre inversion of this formula: /(e) = min[g(A)+Ae]. As g'(X) = —e, 
a simple transformation yields: 
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Figure 2: Surface plot of T curv (a, (3). We took r/7 = 1 for convenience. 




Figure 3: The energy T s i ope associated with the e — > oo divergence. We took r/7 = 1 for 



convenience. 



Figure 4: Exact function g(X) (a = [3 = 1) together with the parabolic Ansatz (dots) based on T c 
(cf. text for details) 




Figure 5: Typical shape of the distribution P(e) for a double well confining potential (See details in 
text), (a): P(e) for different values of r (r = 10, dash-dotted, r = 20, dashed, r = 50, plain), (b): for 
r = 20, comparison between double well potential case (circles) and harmonic potential (crosses) 



